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Advantages of the original symmetrical form of the parametrization of the lepton mixing matrix 
are discussed. It provides a conceptually more transparent description of neutrino oscillations and 
lepton number violating processes like neutrinoless double beta decay, clarifying the significance of 
Dirac and Majorana phases. It is also ideal for parametrizing scenarios with light sterile neutrinos. 
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I. INTRODUCTION 

Since the historic discovery of neutrino oscillations, massive neutrinos currently provide the most direct and testable 
evidence for physics beyond the Standard Model (SM) of particle physics. At low energies, nine (seven) parameters 
must be determined depending on whether neutrinos are Majorana (Dirac) particles Q 1 . Here we tacitly assume the 
former, more general, and theoretically preferred case. Parametrizing the lepton mixing matrix 0, Q in a convenient 
and intuitive manner is very helpful for data handling and interpretation of the physics, such as neutrino oscillation 
searches in upcoming long- baseline experiments 0, Q or searches for neutrinoless double beta decay 0, 0] • While the 
former are sensitive to the Dirac phase, Majorana phases [H, are crucial to describe the latter. 



The Particle Date Group (PDG) has adopted a parametrization of the lepton mixing matrix in which it is a product 



of three consecutive rotations multiplied with a diagonal phase matrix P containing the Majorana phases 12j. The 
mixing matrix can be written as 

U = R 23 (6 23 ;0) R 13 (6 13 ; S) Ri 2 (e 12 ;0) P , (1) 

where Rij(0; ip) is a rotation around the ij'-axis, e.g. 

/ cos 6»i3 sin 0i3 e-" \ 
Rx 3 (0is;6)=\ 10. (2) 
y - sin 0i3 e** cos6>i 3 J 

The position of the Dirac phase S is the convention chosen by the PDG. The two Majorana phases, denoted here a 
and /3, are usually put inside P, to the right of the mixing matrix: P = diag(c M , e 8 ^, 1). However there is no consensus 
notation yet in what concerns the parametrization of these phases, neither for their names (e.g. </>i and 4>2 or ipi and 
ip 2 or a and p, sometimes with a minus sign, sometimes divided by two, etc.), nor for their position within the matrix 
(P = diag(l, c' Q , c i/9 ), or P = diag(e M ,e 1 ^, 1), etc.). The mixing matrix Eq. ([!]) is explicitly given by 



U = UP = 



' ci 2 ci 3 e ia s 12 c 13 e tl3 

-(S12C23 + C12S23S13 e lS ) e M (ci 2 c 23 - si 2 s 2 3Si 3 e lS ) e 10 s 23 c 13 (3) 
\ (S12S23 - ci 2 c 2 3Si 3 e l,5 )e lQ -(ci 2 s 23 + si2C 23 si3e l,5 )e^ 
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1 Many more parameters exist if neutrino masses arise from type-I seesaw, the parametrization in Ref. []J covers all seesaw cases. 
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The elements \U e i\ and \U^\ are known with good accuracy, thus the two "large" angles 812 and 623 are well- 
determined by solar and atmospheric neutrino oscillation data. There is also recent evidence for a nonzero value of 
the "small" element |£7 e 3| from the T2K Collaboration 13[ and the glob al neutrino oscillation data sample 
While there are 8 equivalent ways to parametrize the mixing matrix [16| , the above order of rotation is useful in the 
sense that experimentally the straightforwardly measurable elements \U e \\, |J7 e 3| and {U^l allow to directly extract 
the angles in the three-neutrino lepton mixing matrix 2 . In contrast, the other equivalent parametrizations do not 
share this property. 



II. A MORE CONVENIENT PARAMETRIZATION OF THE MIXING MATRIX 



The above form U is nothing but a re-writing of the "symmetrical" form K proposed in Ref. [l[ (apart from factor 
ordering, which was left unspecified in the original paper). Here we would like to argue in favor of the conceptual 
advantages of the original "symmetrical" presentation of the lepton mixing matrix []J . For the case of three neutrinos 
it is given as: 

K = W 23 (^23 ; 023)^13(6*13; 013)^12(012! 012) , (4) 

where each of the lj's is effectively 2x2, characterized by an angle and a CP phase, e.g. 



/ C13 e 



Wl3 



C13 U e-^ 13 s 13 



10 

\-e^'8i 3 C13 

Explicitly, the symmetrical parametrization of the lepton mixing matrix, K can be written as: 



/ ci 2 ci3 si 2 ci 3 e-^ 12 siae"* 3 



C12C23 - Sl2Sl3S 2 3e 



-i(012 + <fe3-013) Cl3S23e -^23 _ (5) 



\ si2S2 3 e i(012+023) - ci2Si3C 2 3e* 3 -c 12 S2 3 e 1 ^ 3 - s^s^se"^ 12 -^ c 13 c 23 



K 



Here all three CP violating phases are physical [9[ : <pi2 , 023 and tpi 3 . Even though the parametrization is fully 
"symmetric" there is a basic difference between Dirac and Majorana phases. In order to understand this let us use 
the identity [l} 

P^ 1 K P = ^23(6*23; 023 - P) wi 3 (0i3; 013 ~ a) ^12(6*12; 012 + p — a), (6) 
which allows us, up to unphysical phases, to identify the Dirac phase as 

S -H- 013 - 012 - 023 , (7) 

This formula relates the Dirac phase, which denotes the phase responsible for CP violation in neutrino oscillations 
in the PDG parametrization, with the equivalent phase in the symmetrical representation. Note that it obeys field 
rephasing invariance, as it should. Moreover, in contrast to the PDG description, in the symmetrical form CP violation 
in neutrino oscillations is immediately recognized as a three-generation phenomenon involving the phases of all three 
generations 3 , an important conceptual advantage. 



2 See Ref. [17| for a recent application of the other possible parametrizations, and [lgt | for a rare use of the symmetrical parametrization 
we will study in this work. 

3 In this sense the Dirac phase has an intrinsic geometric meaning, like the curl of two vectors or the area of a triangle. 
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Recall that the neutrino oscillation probability for a v a — > up flavor transition is given by 

2 

m 2 

P(y a -> Vp 



= S a p ~ 4J2 » {UriUajUpiUfc} sin 2 

+2^S{[/* i [/ ai ^ i ^.}sin 



Am?- N 

4£ , 



2£ , 



where E is the neutrino energy, L is the distance traveled by the neutrino, and Am 2 



m 2 (m^ being positive 



mass eigenvalues) are the neutrino mass-squared differences. Here 3R and 3 denote real and imaginary parts. For 
three families there is only one independent imaginary part 3 {u^UajUpiUp^, which is responsible for CP violation 
in neutrino oscillations. Comparing this invariant with the PDG and symmetrical parametrizations gives 



j CP - ^{u* el u; 3 u e3 u^} = 

and shows the same result as Eq. ([7]). 



| sin26>i2 sin 26*23 sin26>i 3 cos6»i 3 sin S (PDG) , 

i sin 26*12 sin 26*23 sin26*i3 cos 6*13 sin(</>i3 — ^12 — 023 ) (symmetrical) , 



A. Application to Lepton Number Violating Phenomena 

As well known, there are two conceptually different kinds of CP violating phenomena [l| . In the language of the 
PDG parametrization, one is associated to the "Dirac phase" S and is the exact analogue to the CP phase in the quark 
mixing matrix, responsible for the area of the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle; while the other 
one is associated to the two "Majorana phases" a and /3, which do not show up in neutrino oscillations but do 

affect lepton number violating amplitudes. 

In what follows we will briefly discuss the role of Majoranaphases in determining the rates characterizing neutrinoless 



double beta decay and ncutrino-anti-ncutrino oscillations 

A suitable parametrization of Majorana phases plays a very important role in interpreting the effective mass 
parameter characterizing the amplitude for neutrinoless double beta decay (see Ref. for a recent review). Its 
explicit form reads 



(m) 



U ej m J 



\c 2 2 c 2 3mi e 2 * a + S 2 2 c 2 3 m 2 e 2 ^ + a ? 3 m 3 e 2iS | (PDG) , 
| c i2 c i3 TO i + s i2 c i3 m 2 e 2l< ^ 12 + s 2 3 m 3 e 22 ^ 13 1 (symmetrical) . 



Only the two Majorana phases should appear in (m) [9|. However this is not at all clear in the PDG presentation. 
In contrast, the symmetrical parametrization provides a manifestly transparent description in which only the two 
Majorana phases appear in (m), as it should. Currently nuclear matrix element uncertainties prevent the extraction 
of Majorana phases from neutrinoless double beta decay. However, should these be circumvented and should the 
determination of the Majorana phases become an issue, then the symmetrical parametrization will surely be preferred 
over the PDG one. 

It has long been known that the lepton mixing matrix characterizing the charged current interaction of Majorana 
neutrinos in gauge theories may have complex entries that conserve CP [l9( . These special CP conserving situations are 



associated with Wolfenstein's CP-signs [20(, when neutrino mass states are CP eigenstatcs with CP parity t^cp = it- 
There are four possible inequivalent sign configurations in the sum in Eq. ([S]): (+ + +), (H ), (H h), and (+ + —). 

These are in correspondence to special values of the Majorana phases, namely <j>i2 = cf)i 3 = 0, </>i2 = </>i3 = 7r/2, 
<pi 3 = ir/2 with cf>i2 = 0, and ^13 = with 1^12 = 7r/2, respectively. Majorana phases would also show up in 
processes analogous to neutrinoless double beta decay, such as decays like K + — > tt~ whose amplitude would 

be proportional to ^2 rrii, and have extremely low branching ratios, see Ref. and references therein. 
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Let us now comment on neutrino-anti-neutrino oscillations. A Gedankcncxperiment looking for anti-neutrinos Dp 
in a beam of neutrinos v a has been suggested in Ref. Q in order to clarify the physical nature of Majorana phases at 
the two-generation level. In the three-generation case the probability for such a process is given as 

2 



P(v a -+ up) = i 



03 ' 



-iEit 



1 



Y J U a3 Up 3 U* m U i 



0i ' 



-i{Ej-Ei)t 



(9) 



leading to complicated transition probability expressions, which will not be explicitly given here. The least unrealistic 
channel is represented by v e to D e transitions, because the transition probabilities go with the ratio of mass over 
energy squared, and electron neutrinos can be produced with much lower energy than the other flavors. For three 
families, P(y e — > D e ) depends only on the Majorana phases 0i 2 and 0i 3 for the symmetrical parametrization, whereas 
the PDG case leads to a dependence on all three phases. Note the analogy to the effective mass discussed above. 



B. Application to seesaw and sterile neutrinos 

We now turn to the lepton mixing matrix characterizing models containing gauge singlets such as seesaw models. 
Their most general form was presented within the symmetrical parametrization in Ref. [l| , covering seesaw schemes 
of all types, type-I, type-II, and type-Ill. Since it applies to an arbitrary number m of non-doublet leptons (singlets in 
type-I and II, triplets in type-Ill) this parametrization also covers low-scale in addition to high-scale seesaw schemes. 
To a good approximation, in the standard high-scale seesaw case neutrino oscillations are well-described by the 
simplest unitary form of the lepton mixing matrix used above in Eqs. ([TJ or (j4|). In contrast, for the low-scale seesaw 
schemes 2ll - l24| ne utrino oscillations involve only a unitarity-violating truncation of the full mixing matrix, see, for 



example, Refs. |25|,|26j|. In both cases one has an "effective" neutrino oscillation description with m = 0, i.e. the extra 
neutral states are too heavy to take part in the oscillation phenomena. Since these possibilities have already been 
widely discussed in the literature here we will focus on the alternative possibility that singlets are light enough to 
participate in oscillations, in the s imp le case of m = 1, i.e. one "sterile" neutrino plus three active ones. This possibility 
has recently re-gaii 
We have therefore 



s imp 

0. 



has recently re-gained attention [27j. In terms of the mixing matrix, a useful order of rotation is 34-24-14-23-13-12. 



U = ^34(6*34; 0) ui2i{6%i, 5 2 i) Wi4(0i4; Su) ^23(6*23; 0) u>i 3 (0i3 : 6x3) lo 12 {9 12 \ 0) P , (10) 
with P = diag(e IQ , e 2 ^, e 17 ) in the sense of the PDG description, or in the symmetrical form: 

K = ^34(6*34; 34 ) u; 24 (#24, 024) Wu(014; 01<O ^23(^23; 023) Wi 3 (0 13 ; 013) Wi 2 (012; 012) • (n) 

First note that the number of rotations for 3 + 1 neutrino types (six) is exactly the number of phases (3 Dirac and 3 
Majorana), a characteristic feature of the symmetrical parametrization [l|. 

Consider first the effective mass characterizing neutrinoless double beta decay in this case, which is given as 

(m) - { l Cl2Cl3Cl4 TOl + s i2 c i3 c i4 ™2 c 2lf3 + s? 3 c 2 4 m 3 e 21 ^"- 5 ") + ^ m ^ e -2i6 li j (PDG) , 
~ y \cl 2 cl 3 cj 4 mi + s? 2 c? 3 c? 4 m 2 e 2 ^ 12 + s 2 3 c 2 4 m 3 e 2 ^ 13 + s\ A m 4 e 2 ^ 14 1 (symmetrical) , 

Obviously the symmetrical parametrization has advantages over the one of PDG. Indeed, again, the CP conserving 
cases corresponding to different CP-signs are obtained by choosing the three Majorana phases 0i2,i 3 ,i4 to be 7r/2 or 
zero. 

In what regards oscillations, the three Dirac phases in the PDG description arc <5i 3 , £14 and (524. While there are nine 
[3] independent = S ^U* { U^ - U a j f/gi|, only three independent CP asymmetries P(y a — > vp) — P(yp — > v a ) 
exist. In the symmetrical parametrization, the relevant independent phase combinations appearing in oscillation 
probabilities are 

A23 = 012 + 023 - 013 , 

^124 = 012 + 024 - 014 , (12) 
^134 = 013 + 034 - 014 ■ 
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Each of the phase combinations Iijk, with i < j < k, is seen to span three generations, as necessary for the existence 
of CP violation in neutrino oscillations. Note that there is a fourth possible combination, 1234- However, since 

^123 + Il34 — I\2i = ^234 (13) 



holds, this fourth invariant is not independent. This condition actually implies that I^k is a 2-cocycle |29j. This 
is true for an arbitrary number N s of additional sterile neutrinos: noting that the | N s (N s — 1) rotations between 
sterile neutrinos are unphysical, it is easy to see that for N massive neutrinos, including N s = N — 3 sterile neutrinos, 
there are N — 1 = N s + 2 Majorana phases and 2N — 5 = 2 N s + 1 Dirac phases. Each massless neutrinos results in 
one Majorana phase less. The total number of 3 (N — 2) = 3 (N s + 1) phases corresponds exactly to the number of 
physical rotations, which is \ N (N — 1) — h N s (N s — 1) = 3 (N — 2). The symmetrical parametrization is therefore 
tailor-made also for concisely describing the phenomenology of sterile neutrinos. 

It has been argued that current neutrino data might imply in fact that 2 sterile neutrinos are present f3ol ]. A 
possible order of the nine physical rotations is 



K = U)25 (#25; 025) ^34(^34; 034) W 35 (035; 35 ) W 2 4(#24, 024) 

^23(^23; fe) w 16 (0 15 ; (f>u) u>i 4 (0 14 ; 4>u) wi 3 (0i 3 ; 0i 3 ) ^12(^12; ^12) ■ 
Again, the effective mass characterizing ncutrinoless double beta decay automatically looks straightforward 

/ \ 12222 2 2 2 2 *2i<b-\i ( 2 2 2 2i(Ai^ . 2 2 2i(b-\A • 2 2id>i I 

(m) = c 12 c 13 c 14 c 15 mi + s 12 c 13 c 14 c 15 e 912 + s 13 c u c 15 m 3 e 913 + s 14 c 15 m 4 e <pi4 +s 15 m 5 e <P15 . 



(14) 



involving just the four physical Majorana phases. Regarding CP violation in oscillations, there are 
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10 — 3 = 7 different I^k combinations with i < j < k and i, j € {1, 2, 3, 4, 5}, where the subtraction of 3 stems from 
the cases which have both 4 and 5 in ijk: 




-fl23j ^124, Il25, ^134j ^135, ^234, ^235 ■ (16) 

There also exist I 1—3 = 5 — 3 = 2 "sumrulcs" in analogy to Eq. (Ti"3| , namely 

^123 + -^134 — I\2i = ^234 and /123 + -/135 — I125 = ^235 • (17) 

Hence, at the end there are 5 physical Dirac CP violating phase combinations affecting neutrino oscillation probabil- 
ities, for instance one could choose 



I\23 = < 


j>X2 + 


023 - 


013 , 


(18) 


I\2i — < 


i>u + 


024 - 


014 , 


(19) 


1 134 = 1 


ha + 


034 - 


014 , 


(20) 


I\25 = 1 


i>i2 + 


025 ~ 


015 , 


(21) 


^135 = 


013 ■+ 


" 035 - 


015 • 


(22) 



The phase relevant for CP violation in the short-baseline oscillations v e O-i^ sector is 3 {U* 4 U* 5 J/es} <x sin(0i 4 - 

015 - 024 + 025) = Shl(/i25 - /l2 4 ). 

The generalization to N s sterile neutrinos is now clear: In more mathematical language, adopted from Rcf. 
one may define the operator 6 (with S 2 = 0) such that (cf. Eq. (TT31) ') 

5/1234 = F% = /i2 3 + J 134 - /l24 - ^234 = . (23) 



For general i,j,k,l with i < j < k < I one has Sliju = 0. For five active generations one has ten I^j- and five 

pressed by the other four, for instai 

P (4) _ p(4) p (4) p (4) p (4) , p _„ (0A , 

^2345 - ^1235 + ^1345 ~ -^1234 _ ^1245 0r <^12345 - U . (24 J 



Fijkt = 0, one of which can be expressed by the other four, for instance 
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Therefore, the standard result of 10 — (5 — 1) = 6 Dirac phases is obtained. If all generations were active, then the 
number of independent phase combinations is 







\(N-l)(N 



(25) 



One simply subtracts the number of | N s (N s — 1) = | (JV — 3) (N — 4) unphysical rotations from this result to obtain 
the already quoted 2N — 5 = 2N a + l Dirac phases. The first binomial in Eq. (|25l) is the number of ijk combinations 
with i < j < k, the second binomial is the number of sumrules between them, and the third binomial describes 

f:,\ 

the linear relations existing between the sumrules. For instance, 5 active generations have 



combinations with i < j < k and i, j G {1, 2, 3, 4, 5}. A number of 
relation between them. 



5 sumrules exist, which 



10 different i^-fc 
^|=1 linear 



III. FINAL REMARKS 

In summary, the issue of a proper parametrization scheme for the lepton mixing matrix will become relevant as 
experiments reach sensitivity to CP violation, either of Dirac or Majorana type. Even more so, if present indications for 
sterile neutrinos are confirmed in upcoming experiments. While the form given in the PDG is just a re-writing of the 
symmetrical form proposed in Ref. |l| , we have advocated here the conceptual advantages of the original symmetrical 
presentation for the description of neutrino oscillations and, especially, of lepton number violating processes like 
neutrinoless double beta decay. 
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